Multi-scale analysis of the electromagnetic self-force in a weak gravitational field 
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We examine the motion of a charged particle in a weak gravitational field. In addition to the 
Newtonian gravity exerted by a large central body, the particle is subjected to an electromagnetic 
self-force that contains both a conservative piece and a radiation-reaction piece. This toy problem 
shares many of the features of the strong-field gravitational self-force problem, and it is sufficiently 
simple that it can be solved exactly with numerical methods, and approximately with analytical 
methods. We submit the equations of motion to a multi-scale analysis, and we examine the roles of 
the conservative and radiation-reaction pieces of the self-force. We show that the radiation-reaction 
force drives secular changes in the orbit's semilatus rectum and eccentricity, while the conservative 
force drives a secular regression of the periapsis and affects the orbital time function; neglect of the 
conservative term can hence give rise to an important phasing error. We next examine what might 
be required in the formulation of a reliable secular approximation for the orbital evolution; this 
would capture all secular changes in the orbit and discard all irrelevant oscillations. We conclude 
that such an approximation would be very difficult to formulate without prior knowledge of the 
exact solution. 



I. INTRODUCTION 



The gravitational inspiral of a solar-mass compact ob- 
ject into a massive black hole residing in a galactic center 
has been identified as one of the most promising sources 
of gravitational waves for the Laser Interferometer Space 
Antenna [H . The need for accurate theoretical models of 
the expected signal, for the purposes of signal detection 
and source identification, has motivated an intense effort 
from many workers to determine the motion of the small 
body in the field of the large black hole. This is done 
in a treatment that goes beyond the geodesic approxi- 
mation and takes into account the body's own gravita- 
tional field, which is a small perturbation over the field 
of the black hole. In this treatment the small body can 
be described as moving on an accelerated world line in 
the background spacetime of the large black hole; the 
body is said to move under the influence of its own grav- 
itational self- force 0, and this force is derived from 
the retarded gravitational perturbation produced by the 
moving body. For a review of the self-force formalism, see 
Ref. Q and the special issue of Classical and Quantum 
Gravity devoted to this topic fs*] . 

The concrete evaluation of the gravitational self-force 
acting on a small body moving in the Kerr spacetime is 
a challenging project that has not yet been completed 
(although progress has been steady). Given the severity 
of the chall enge, a number of authors [1, 0, H, d, 
[13, [Hj Hj llSl . [l6| have attempted to formulate various 
simple schemes that would allow them to reproduce the 
effects of the self-force; their hope is that these schemes 
will be simple enough for rapid implementation in nu- 
merical codes, and accurate enough to describe faithfully 
the orbital evolution of a body subjected to its gravi- 
tational self-force. One such scheme is Mino 's radiative 
approximation 0, 0j @|, which is based on an approxi- 
mate self-force constructed from the half-retarded minus 
half-advanced gravitational perturbation associated with 



the moving body. Mino proved that while his version of 
the self-force neglects all conservative corrections to the 
motion, it correctly accounts for the long-term dissipa- 
tive effects associated with the true self-force. This led to 
the widespread belief that all long-term secular changes 
in the orbital motion would be captured by the radiative 
approximation, and that conservative effects would pro- 
duce only short-term changes that would not accumulate 
in the long run. The simplicity of Mino's scheme made 
it attractive, and it was adopted in a number of works 
that aimed to model the inspiral of a small body into a 
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rapidly rotating black hole 

Mino's radiative approximation was criticized, how- 
ever, in an earlier work by the present authors [T7|, here- 
after referred to as "paper I." Building on an analogy be- 
tween the gravitational self-force and its electromagnetic 
counterpart, we showed that conservative terms in the 
true self-force do lead to long-term secular changes in the 
orbital motion. These changes are not captured by the 
radiative approximation, and we concluded that Mino's 
scheme has severe limitations. This conclusion was sup- 
ported by a recent analysis by Drasco and Hughes 
and the general attitude currently is that while the radia- 
tive approximation may be useful to generate template 
waveforms for signal detection, it is probably insufficient 
for reliable parameter estimation. 

Our purpose in this paper is to revisit the analysis pre- 
sented in paper I. There are three reasons for reopening 
the case. The first is that our original analysis of the elec- 
tromagnetic self-force employed rather crude mathemat- 
ical tools, and we wish to present here a more thorough 
and rigorous treatment. The second is that the main 
source of discrepancy between the effects of the radiative 
self- force and those of the true self- force was not correctly 
identified in paper I. In the original paper we claim that 
the discrepancy is mainly due to the secular regression 
in the orbit's periapsis, an effect that is produced by the 
true self-force but not accounted for by the radiative ap- 
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proximation. In this paper we show that while this is 
indeed a source of discrepancy, it is not the most impor- 
tant one. As we shall explain in Sec. IV C, the most 
important conservative effect is actually associated with 
the time function on the orbit. The third reason is that 
we wish to introduce here a clear distinction between the 
radiative approximation to the self-force and the notion 
of a secular approximation to an orbital evolution; our 
secular approximation is a specific implementation of the 
general idea of capturing the long-term orbital evolution 
through an adiabatic approximation that allows the orbit 
to evolve slowly. The phrases "radiative approximation" 
and "adiabatic approximation" are used synonymously 
in paper I (and indeed, in most of the literature on this 
topic), but we feel that this is a highly misleading prac- 
tice. A large portion of this paper is devoted to the task 
of identifying what should be required of a good secu- 
lar approximation, and we shall see that the radiative 
approximation does not meet those requirements. 

The precise meaning of an adiabatic approximation is 
somewhat ambiguous in the literature. In all cases, the 
basic assumption is that the secular effects of the self- 
force occur on a time scale that is long compared to the 
orbital period. From this assumption, numerous approx- 
imations have been formulated: (1) Since the particle's 
orbit deviates only slowly from geodesic motion, the self- 
force can be calculated as if the particle travels on a 
geodesic (or, in the post-Newtonian case, the radiation 
reaction can be calculated as if the particle's dynamics 
were conservative); (2) since the radiation-reaction time 
scale is much longer than the orbital period, periodic 
effects can be neglected; and (3), based on various argu- 
ments, conservative effects can be neglected. Although 
each one of these three approximations has been called 
an adiabatic approximation, we believe that they should 
be distinguished from one another. To discuss the first 
approximation is beyond the scope of this paper. We 
focus instead on the latter two approximations: number 
(2) above, which we call the secular approximation, and 
which neglects periodic effects; and number (3), which we 
shall call the radiative approximation, and which neglects 
conservative effects. 

The main idea behind the construction of a secular 
approximation is the following. We consider an orbital 
evolution under the action of a self-force, and we wish to 
simplify the equations of motion in such a way that the 
long-term, secular changes will be captured, at the cost 
of discarding irrelevant, short-term effects. Suppose that 
we describe the orbital evolution in terms of a set of or- 
bital elements I'^it), where A is an index that labels each 
element. (This description is introduced in Sec. Ill, and 
explained fully in Appendices A and B.) The orbital ele- 
ments would be constant in the absence of a perturbing 
force, but they evolve in time as a result of the force's ac- 
tion. It is expected that each orbital element will display 
a behavior that can be decomposed into a secular change 
that accumulates monotonically over time, and an oscilla- 
tion that averages to zero in the long run. We thus write 



(t) = (t) -I- /^(, (t) , and a secular approximation for 
the orbital elements would keep the secular terms and 
discard the oscillations. We would write I^{t) ~ 4'cc(0i 
and seek a method to obtain I^^i^) in the most direct 
and convenient way possible. Presuming that this must 
be done in a context in which the exact solution I^{t) 
would be too difficult to obtain, we would seek to for- 
mulate equations of motion directly for I^^i^), and we 
would hope that those equations are sufficiently simple 
that a solution could easily be found (analytically or nu- 
merically) . This is the main idea, and the task ahead ap- 
pears to be clearly identified. But to turn the idea into 
a precise algorithm may not be easy. To illustrate the 
difficulty we shall examine, in a specific context in which 
we can make progress (the electromagnetic self-force of 
paper I), what would be required in the construction of 
a secular approximation for the orbital evolution. 

The secular approximation is logically distinct from 
the radiative approximation, in which the true self-force 
is truncated so as to discard all conservative terms. In 
the radiative approximation, one writes I"^{t) ~ ^^dC^)' 
and one calculates I^^{t) on the basis of the truncated 
self- force. It is known, as Mino has shown [1, 0, [1| , that 
the radiative self-force correctly accounts for the long- 
term, dissipative changes in the orbital elements. If it 
correctly produced the long-term, conservative changes 
as well, we would conclude that the radiative approxima- 
tion captures the idea of a secular approximation. But, 
as we have shown in paper I [17], and as we intend to 
show even more convincingly here, the radiative approx- 
imation fails to account for secular changes in I"^{t) that 
are produced by the conservative piece of the self-force. 
The radiative and secular approximations are therefore 
distinct, and we consider it important to distinguish these 
terms carefully. 

We will examine the limitations of the radiative ap- 
proximation, and attempt to construct a faithful secular 
approximation, in the specific context of an electromag- 
netic self-force acting on a charged particle moving in a 
weak gravitational field. The motion of the particle is 
governed by the equations 

a = 9 + /self, (1-1) 

where a = d^r/dt^ is the particle's acceleration vector, 
g — ~Mr/r^ is the Newtonian gravitational field of a 
body of mass M , and 

q^M 2q^dg 

/self = Ac -I- A„.- — — (1.2) 

/i r"' 6 ^ at 

is the electromagnetic self-force divided by the particle's 
mass /i. (The particle is treated as a test mass. In this 
treatment, the central body does not move, its mass M 
is the system's total mass, and /z is the s yst em's reduced 
mass.) This self-force was calculated [Tsl . ITgj for the 
weakly curved spacetime produced by the central body, 
assuming that the charged particle moves slowly. Here 
q is the particle's charge, and r{t) is its position vector 
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relative to the central body; we have also introduced the 
distance r = \r\ and the unit vector r — r/r. The con- 
stants Ac and Air in Eq. (|1.2p are both equal to unity; 
they serve to remind us that the first term in Eq. (|1.2p 
is the conservative piece of the self-force, while the sec- 
ond term is the dissipative (or radiation-reaction) piece. 
By keeping these constants in our calculations we will be 
able to distinguish conservative effects from dissipative 
effects; the radiative approximation is obtained by set- 
ting Ac = and keeping X„ = 1. Throughout the paper 
we use the usual vectorial notation of three-dimensional 
flat space, and we work in units such that G = c = 1. 

We work in the specific context of Eqs. and (|1.2[) 
for two reasons. First, our toy problem is an actual ex- 
ample of a self-force acting on an orbiting body. While 
the force has an electromagnetic origin instead of a grav- 
itational origin, and while the motion takes place in a 
weak (Newtonian) gravitational field instead of a strong 
field, the self-force of Eq. (I1.2|l nevertheless contains con- 
servative and dissipative terms that will have different 
effects on the orbital motion. In the usual Lorenz gauge, 
the gravitational self-force in strong fields will also con- 
tain conservative and dissipative pieces. In addition, each 
self-force comes with a similar post-Newtonian counting. 
From Eq. (jl.2p we see that the conservative term in the 
electromagnetic self- force is a correction of order q"^ /{pr) 
relative to g, and taking q to be of order fi, we recognize 
this as a correction of first post-Newtonian (Ipn) order; 
the dissipative term is a correction of order q^v/{^r), 
where v is the orbital velocity, and we recognize this as a 
correction of 1.5pn order. On the other hand, in a post- 
Newtonian context the gravitational self-force presents 
conservative pieces at orders Opn, Ipn, 2pn, 3pn, and 
so on, and dissipative pieces at orders 2.5pn, 3.5pn, and 
so on. While the post-Newtonian counting is not iden- 
tical, in each case we have dominance of the conserva- 
tive effects over the dissipative effects, and all in all, this 
gives us good reasons to believe that the electromagnetic 
problem captures the essential physics of the more com- 
plicated, gravitational problem when it is formulated in 
the usual Lorenz gauge. 

Second, Eq. is far simpler than the realistic self- 
force equation (for which the force must be obtained nu- 
merically), and this permits a very thorough and rigor- 
ous mathematical analysis. We shall therefore be able 
to extract very precise consequences of Eqs. (jl.ip and 
(II. 2p . and examine closely the issues that concern us re- 
garding the radiative and secular approximations. The 
simple mathematics of the toy problem will allow us to 
draw firm and clear conclusions, and the proximity of its 
physics to that of the realistic problem will give us confi- 
dence that these conclusions extend from the toy problem 
to the realistic case. 

We begin in Sec. II with a simple illustration of the 
themes to be explored in this paper. The mathematical 
analysis of Eqs. (|l.ip and (|1.2p is carried out in Sec. Ill, in 
the framework of osculating orbital elements summarized 
in Appendices A and B. The mathematical details are 



presented in Sec. Ill with minimum commentary, but we 
present a full discussion of our results in Sec. IV. We 
summarize our conclusions in Sec. V. 



II. ILLUSTRATION 

Before we proceed with our mathematical analysis of 
Eqs. (|l.ip and (|1.2p . it is helpful to describe a very simple 
problem that illustrates rather well the issues we shall 
encounter. 

Suppose that we are interested in a quantity q{t) that 
is governed by the system of dynamical equations 

dq dt 

— =ei-e2SmA, — = H- cos A -I- €3, (2.1) 
aX d\ 

where A is a running parameter, and ei, 62, and 63 are 
small constants. The differential equations come with the 
initial conditions 

g(0) = 1, m = 0. (2.2) 

In this example, q{t) is analogous to the set of orbital 
elements I^{t) that were introduced previously, and the 
equation for dq/dX is analogous to Eq. p.lOp below, with 
ei and €2 playing the roles of e„ and ec, respectively. The 
quantity q is constant when Ci = £2 = (the unperturbed 
situation), and it acquires a time dependence when the 
perturbation is turned on. The equation for dt/dX is 
analogous to the first line of Eq. (|3.34p below, with 
playing the role of ec. As we explain in Sec. Ill C, the 
first term proportional to Cc in Eq. p.34p is generated by 
oscillatory terms in the orbital elements and oscillatory 
terms in the unperturbed equation for t; these combine 
to give rise to a secular term in the perturbed equation. 
To ignore the oscillations would produce the significant 
mistake of dropping the £3 term in Eq. (|2.ip . The vari- 
able A gives a convenient parameterization of the motion; 
its use is motivated by the fact that the "force" can be 
expressed as a simple function of A, while it would be 
very difficult to express it in terms of the time variable t. 
The parameter A is analogous to the orbital parameter (jj 
that will be introduced in Sec. III. (It is also analogous 
to "Mino time" [l, [l^l , a convenient parameterization of 
geodesies in Kerr spacetime.) 

The exact solution to the system of equations is 

q{X) = l + eiA + e2(cosA-l), (2.3) 
<(A) = A + sinA + e3A. (2.4) 

We see that ei produces a secular growth in q{X), 62 is 
associated with oscillations, and £3 produces a secular 
drift in the time function t{X). 

Suppose that we are interested only in the long-term, 
secular changes in q, and that we wish to construct a 
secular approximation for it. Because we have access to 
the exact solution, it is a simple matter to remove the 



4 



oscillations by subjecting it to an averaging procedure. 
Our first option is to introduce 

{q)x:=^ r\{X')d\' (2.5) 

and to define the secular approximation as gsoc = ('?)a- 
Our choice here is therefore to remove the oscillations 
with respect to A, and a calculation based on Eqs. (|2.3p 
and (|2.5|) gives 

(<z)a = l + eiA-e2 + 0(e'). (2.6) 

This version of the secular approximation is a solution to 
the modified differential equation 

^{q)x = e, + 0{^) (2.7) 
with the modified initial condition 

(g)A(0) = l-e2 + O(e'). (2.8) 

If we did not have access to the exact solution, we might 
still have guessed that the correct differential equation 
for qscc is Eq. (|2.7p . because it can be obtained directly 
from Eq. (|2.ip by averaging over the oscillatory term. 
But we would be hard pressed to guess that the correct 
initial condition is given by Eq. (|2.8p . Using the approxi- 
mate differential equation with the exact initial condition 
Qscc(O) = 1 would produce a function that is offset by 62 
relative to {q)\- 

Our first message is that a faithful secular approxi- 
mation can be based on an averaged version of the dif- 
ferential equation, but that it must come also with a 
corresponding change of initial condition. To obtain the 
approximate differential equation might be easy, but to 
identify the correct initial condition is impossible when 
the exact solution is unknown. The formulation of a sec- 
ular approximation therefore suffers from an ambiguity 
regarding the correct choice of initial condition. In this 
example the consequence of missing the £2 term in the 
initial condition is not severe: The difference between 
the solutions 1 -I- eiA and 1 -I- eiA — £2 becomes relatively 
small as A increases and each solution grows secularly. In 
other situations, however, the difference in initial condi- 
tions could lead to more serious discrepancies. 

In Eq. (|2.5p we removed the oscillations of the exact 
solution by averaging over the parameter A. Because the 
observer might be more interested in the time-behavior 
of the function g, an alternative choice is to perform the 
averaging over t instead of A. And since t{X) contains 
oscillations, it should be expected that this alternative 
method of averaging will lead to a distinct formulation 
of the secular approximation. Our second option is there- 
fore to introduce 

iii:<,{>.'midy)dx 



and to define version 2 of the secular approximation as 
9soc = {<i)t- A calculation based on Eqs. (|2.ip . (|2.3p . and 
gives 

{q)t = 1 + ei(A - sin A) - ^£3 + 0{e^). (2.10) 
This is a solution to the modified differential equation 

-^(g), = + cos A) + 0(£2) (2.11) 
dA 

and the modified initial condition 

{q)m^l-\e2 + 0{e^). (2.12) 

Here the situation is more interesting. If we did not have 
access to the exact solution, we would never have guessed 
that the correct differential equation for the secular ap- 
proximation is Eq. (|2.1ip . and we would also never have 
arrived at Eq. ([^1^ . 

Our second message is that this new secular approx- 
imation (version 2, which removes the oscillations in t 
instead of the oscillations in A) must be based on an 
approximate differential equation and an approximate 
initial condition that are impossible to identify without 
knowing the solution to the exact problem. The ambigu- 
ity of the first method extends from the choice of initial 
condition to the specification of the differential equation. 

Our third message is that while the idea of formulat- 
ing a secular approximation is clear enough, it is difficult 
to turn it into a precise algorithm. To remove the os- 
cillations of an exact solution is easy enough. But to 
reformulate the system of differential equations and ini- 
tial conditions into a set of approximate equations that 
would achieve the same result is difhcult; it might well 
be impossible in most cases. 

Our fourth message is concerned with the ana- 
logue here of formulating a radiative approximation to 
Eqs. (|2.ip . This is obtained by setting £2 = £3 = while 
leaving £1 unchanged. This produces the functions 

grad(A) = l + £iA, irad(A) = A + sinA. (2.13) 

After a long time, when A 3> 1, 5rad(A) becomes very 
nearly equal to q{\), and the radiative approximation is 
accurate when q is expressed in terms of the orbital pa- 
rameter. At late times, however, we have that trad (A) — A 
while t(A) ~ (1 -|- £3)A, and we see that the radiative ap- 
proximation produces a shift in the time function that 
becomes important when A increases beyond l/£3; we 
shall see that this feature is present also in the context of 
the electromagnetic self-force in Sec. III. When q is ex- 
pressed as a function of time, we get that (j'rad(t) — l+ei^, 
while the exact solution behaves as q{t) ~ H-£it/(l + £3). 
The difference is equal to £i£3i/(l -I- £3). While this ap- 
pears to be small because of the first factor of order £^ , it 
is steadily growing because of the additional factor of t. 
The radiative approximation, therefore, produces a sec- 
ular drift in the time function, and a corresponding drift 
in q. 
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III. MULTI-SCALE ANALYSIS OF THE 
ELECTROMAGNETIC SELF-FORCE 

We now proceed with our mathematical analysis of 
Eqs. (jl.ip and (|1.2p . This section, unlike all others in 
this paper, is highly technical, and we intend to deal with 
the technical issues while keeping the commentary to a 
minimum. The implications of our results, in the light 
of the themes introduced in Sec. I, will be fully detailed 
in Sec. IV. The reader who may not wish to delve into 
the technical details, and who would prefer to pick up 
the story where we left off at the end of Sec. II, can omit 
reading this section and proceed directly to Sec. IV. 



A. System of equations 

We wish to integrate the equations of motion (|l.ip for 
the electromagnetic self- force of Eq. (|1.2p . We shall do so 
by employing the method of osculating orbital elements 
developed in Appendix B. 

The starting point of the method is the unperturbed 
situation described by the equations a = g, in which 
the particle follows a Keplerian orbit characterized by 
a number of orbital elements. (Kepler's problem is re- 
viewed in Appendix A.) The orbital elements are con- 
stants of the Keplerian motion; they are related to the 
initial conditions placed on the particle's position and ve- 
locity vectors, but they are defined so as to provide the 
most useful information regarding the geometric proper- 
ties of the orbit. The elliptical shape of the Keplerian 
orbit is described by 



be described by Eqs. (|3.ip and (|3.2p . but the orbital ele- 
ments (p, e,w) acquire a 0-dependence that accounts for 
the perturbation. Their evolution equations are given by 
Eqs. (|B16[) - (|B18[) in Appendix B. They rely on a decom- 
position of the self-force according to /self — Rr + S(f), 
with R denoting its radial component and S its tangen- 
tial component; the unit vectors r and point in the 
directions of increasing r and 0, respectively. The self- 
force does not contain a component normal to the orbital 
plane, and indeed, our version of the method of osculat- 
ing elements is restricted to perturbing forces that are 
tangent to the plane. 

From the expression given in Eq. (|1.2p . we find that the 
radial and tangential components of the self-force are 



„ g^M/ 4^ A 1 
R=- I Ac + -A„r- -r 



and 



(3.3) 



(3.4) 



respectively. Here an overdot indicates differentiation 
with respect to t, and time derivatives can be converted 
into (/)-derivatives by involving Eq. (j3.2p . After differen- 
tiating Eq. (|3.ip to obtain r', we find that Eqs. (|3.3p and 
become 



R 



q^M (1 4- ecf 



4^ M 
Ac + o AitW — es 
3 \ P 



(3.5) 



and 



r(0) 



1 -I- ecos(0 — uj) 



(3.1) 



where r is the distance between the particle and the cen- 
tral mass M, and is the longitude. The orbital ele- 
ments are the semilatus rectum p, the eccentricity e, and 
the longitude at periapsis u. The elements p and e deter- 
mine on which ellipse the particle is moving, and we shall 
call them the principal orbital elements. The element uj 
determines the particle's initial position on the selected 
ellipse, and we shall refer to it as a positional orbital el- 
ement. The position of the particle as a function of time 
is determined by integrating 
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M [1 + ecos{4> - Lu)Y 



(3.2) 



for the time function t{(f>); the prime indicates differenti- 
ation with respect to (/>. The motion in the orbital plane 
is then fully described (in parametric form) by the func- 
tions r{(p) and t{(p). It is an important fact that Eq. p.2p 
does not admit a closed-form solution; a convenient way 
to handle it is by straightforward numerical integration. 

We next move on to the equations a = g + fscu and 
a description of the perturbed motion. In the method 
of osculating orbital elements, the motion continues to 



S ■ 



q^M (l + ec)"^ 




(3.6) 



where c := cos(0 — w) and s := sin((/) — w). These expres- 
sions are ready to be inserted within Eqs. (|B16p - (|B18p . 

The evolution equations come with the initial condi- 
tions 



P(0 = 0) 
e(0 = 0) 

u;{(j) = 0) 
i(0 = 0) 



P , 
e*, 

u* = 0, 
r = 0; 



(3.7) 



the values selected for uj* and t* produce no loss of gen- 
erality. To facilitate the integrations we introduce the 
dimensionless semilatus rectum p and dimensionless time 
t, as well as the dimensionless parameters Ec and e^. that 
characterize the strength of the perturbing force. These 
are defined by 
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possible, and for this purpose we shall construct approx- 
imate solutions to these equations, taking advantage of 
the fact that and are small. This will be carried 
out in the following subsections. The approximation we 
shall construct is distinct from the secular and radiative 
approximations considered in Sec. I; we shall refer to it 
as the multi-scale approximation. We shall demonstrate 
that our multi-scale approximation faithfully reproduces 
the numerical results at all points on the orbit, up to 
a time at which terms of second order in Cc and e„ be- 
come important. From the multi-scale approximation we 
shall be able to construct secular and radiative approxi- 
mations, and we shall be able to ascertain their accuracy. 



FIG. 1: Orbital evolution under the action of the electro- 
magnetic self-force. We set M — 1, q^/fJ, = 0.05, and the 
orbital elements are integrated from the initial conditions 
p* = 20, e* = 0.8, and ui* — 0. With these choices we 
have ec = 2.500 x 10"^ and e„. = 3.727 x 10"*. The motion 
proceeds in the counterclockwise direction and is followed for 
50 orbital cycles {0 < (j) < IOOtt). A the end of the integration 
ec<j) ~ 0.7854. The orbit is displayed in the x-y plane, with 
X = r cos (j) and y — r sin (f>. In the course of its evolution 
the orbit becomes smaller, more circular, and its periapsis 
regresses. 




The final form of the evolution equations is 
1 + ec 



,1/2 ' 



(3.8) 



(3.9) 



s{l + ec) , (l + ec)(e-2c-3ec^) 
e = Cc he„ — , (3.10) 

c(l + ec) (1 + ec)s{2 + Sec) 



ep 



ep 



3/2 



^3/2 



(1 + ecy 



(3.12) 



where c = cos(0 — uj) and s — sm{(/) — to). Integration 
proceeds from the initial values p(0 = 0) = 1, e(0 = 0) = 
e*, a;((/) = 0) = 0, and t(0 = 0) = 0. We shall assume 
that Cc and e„ are small throughout the evolution. In 
spite of the fact that err is smaller than Cc by a factor of 
order y/M/p* <C 1, we shall formally treat them as being 
of the same order of magnitude. 

Equations p.9p - p.l2p can easily be integrated numer- 
ically, and the orbital motion reconstructed by insert- 
ing the solutions within Eq. (|3.1[) . The result of such a 
numerical integration is presented in Fig. 1. Our goal, 
however, is to obtain as much analytical information as 



B. Multi-scale approximation: orbital elements 

The action of the electromagnetic self-force causes the 
orbital elements p, e, and uj to acquire a ^-dependence 
governed by Eqs. (|3.9p - (|3.1ip . The corrections to the 
unperturbed solutions p = 1, e = e*, and to — can 
be separated into two classes: secular terms that grow 
monotonically with and nonsecular terms that oscil- 
late and average to zero over a complete orbital cycle 
(0 < < 2ti). To capture these different behaviors, we 
need an approximation method that has the capability 
of producing a solution that stays accurate over a long 
interval < < (/)rnax, with 0rnax of the order of e~^, 
where e is the overall smallness parameter of the prob- 
lem; and in this interval, the difference between the exact 
and approximate solutions must be uniformly of order 
e^. (Because we have introduced two such parameters, 
we shall write Ec = Cce, €„■ = CrrC and consider Cc, Crr 
to be quantities of order unity.) These requirements rule 
out a simple-minded expansion in powers of e, because 
this method would give rise to a solution that is accurate 
only for £0^1. We adopt instead a multi-scale analysis 
(see, for example. Chapter 11 of Ref. j2]|). 

In a multi-scale expansion one introduces a dependence 
on a "long-scale" variable z := e0 in addition to the 
dependence on the "short-scale" variable 0. We write 



P 

e 

UJ 



Po(2;) + epi(z, ( 
eo(z) -I- eei(z,(; 
uJo{z) + euji{z, 



(3.13) 
(3.14) 
(3.15) 



and we seek to isolate all secular changes within the 
zeroth-order quantities, and to make all first-order quan- 
tities purely oscillatory. We use the chain rule 

d<b dz 



to evaluate the total derivative with respect to 
function /(z, cj)). 



of a 
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To proceed we substitute Eqs. (|3.13p - (|3.15p into Eqs. p.9p - l|3.1ip and obtain, to first order in e, 

2e„. 



dz dcj) 

deo dei 
dz 86 



dujQ duji 
dz d(j) 



1/2 

Po 



1 + eo cos v] , 



Po 



1 



sinw +2^0 sin2w 



3/2 

Po 

'2po 

3/2 

Po 



3 1 5 3 ' 

-eo + ^(8 + Seq) cos w + -eo cos2u + -Eq cos3w 



1 H cos t; + cos 2v 

eo 

8 + 3e^ . 5 . „ 3 . ^ ■ 

sm ti H — sm 2w H — eg sm 3v 

4eo 2 4 



(3.16) 



(3.17) 



(3.18) 



where v := (j) — loq. It is easy to recognize the terms on the right-hand sides that drive the secular changes in the 
orbital elements. We isolate these changes by setting 



dpo 2e,. 
dz 



1/2 ' 

Po 



deo SerrCo 



dujQ 
dz 



Be 

'2po' 



The nonsecular (oscillatory) corrections are then obtained by integrating 

26^60 



dei _ Cc 
dcj) Po 



1/2 

Po 



■ COS?;, 



1 



sinw +2^0 sin 2v 



1 5 3' 

- (8 + Sep) cos V + -eo cos 2v + -Cq cos 3v 



duJi 



3/2 

Po 

ec 
2po 

3/2 

Po 



— COS V + COS 2t; 
eo 

8 + 3e;" 



sin u H — sin 2v -\ — eo sin 3v 
4eo 2 4 



We must impose the initial conditions po + epi + ■ • • = 1, eo + eei + ■■■ = €*, and loq + emi + ■ 



(3.19) 
(3.20) 
(3.21) 

(3.22) 

when 6 = 0. 



The general solutions to Eqs. (|3.19p are 

Po = a(l - 3err^/a^^^)^^^, 
eo = b{l~3e,,z/a^/^f\ 



(3.23) 
(3.24) 



UJQ 



ecfl 



1/2 . 



2e 



1- (l-3e„^/a=^/')'^'^l, (3.25) 



where a, b, and c are constants of integration that will 
be determined. 

The purely oscillatory solutions to Eqs. ([3?20l) - ([3?22l) 
are 



ePi 



eei = 



2ei.reo . 
Po 



(3.26) 



Po 



1 

cosTj + ^eo cos 2?; 



A 3/2 

4Po 



'< + 5eo) sin v + 5eo sin 2v 



etoi = 

Po 



eo sin 3t; 



1 . 1 . 

— sin V -\ — sin 2v 

eo 4 



(3.27) 



A 3/2 

4Po 



8 + 3eg 



eo 

Co COS 3v 



cos V + 5 cos 2v 



(3.28) 



We recall that v = (f> — luq- 

To relate the constants a, b, and c to the initial con- 
ditions we note that Po(0) = a and epi(O) = 0; we 
therefore have a = p* = 1. Similarly, we note that 
eo(0) = b and eei(O) = — iec(4 + b); we therefore have 
b = e* -I- iec(4 + e*) + 0(e2). Finally, we note that 
wo(0) = c and ewi(O) = err(8 + 5b + 462)/(46); we there- 
fore have c = -e„(8 -t- 5e* 4e*2)/(4e*) ©(e^). Making 
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FIG. 2: Multi-scale approximation for p{(p/2n) compared with 
exact numerical results. The numerical conditions are the 
same as in Fig. 1. The solid curve in red shows the exact 
evolution as computed numerically. The dotted curve in blue 
shows the evolution as predicted by the multi-scale approxi- 
mation, which includes a secular term as well as oscillations. 
The dashed curve in green shows the secular piece of the 
multi-scale approximation. The large panel shows the entire 
evolution from (fi — to (fi — IOOtt. The first inset (bottom 
left) shows the evolution in the small interval 6 < (f)/{2n) < 9; 
early in the evolution the multi-scale approximation is ex- 
tremely accurate. The second inset (top right) shows the 
evolution in the small interval 41 < (/)/(27r) < 44; here the 
multi-scale approximation is less accurate, because e0 has be- 
come comparable to unity. 



these substitutions in Eqs. p.22p - (|3.24p gives 

,2/3 



Po 

eo 



= l-3e. 



1 



l-3e. 



.1/2 



Wo = -Cr 



2e„ 



4e* 
5e* -t-4e*2 
4e* 

1 - (l - Serr-} 



,1/3- 



(3.29) 
(3.30) 

(3.31) 



These expressions describe the secular changes in the or- 
bital elements. Equations p.26p - p.28p . on the other 
hand, describe the nonsecular (oscillatory) changes. All 
together, these results give us the desired multi-scale ap- 
proximation for the orbital elements.^ We compare the 



FIG. 3: Multi-scale approximation for e(0/27r) compared with 
exact numerical results. The numerical conditions are the 
same as in Fig. 1, and the caption of Fig. 2 provides the 
relevant details. 




20 30 
(p / (2 Jl) 

FIG. 4: Multi-scale approximation for lu[4>/'2h) compared 
with exact numerical results. The numerical conditions are 
the same as in Fig. 1, and the caption of Fig. 2 provides the 
relevant details. 



approximations with exact numerical results in Figs. 2, 
3, and 4. 



We take this opportunity to make an observation. We notice 
from Eq. 113.291 1 that when is comparable to unity, the dissi- 
pative term in the self-force produces a change in p that is also 
of order unity. In this calculation the radiation-reaction force 
is linear in e, and there are no corrections of order . If such 
corrections were present, however, they would produce an addi- 
tional change of order e in p. Next we notice from Eq. I I3.3UI I 
that the conservative term in the self-force produces a change of 
order e in e, in addition to the change of order unity that comes 
from the radiation-reaction force. We conclude that second-order 



terms in the radiation-reaction force would produce effects that 
scale with the same power of e as those produced by the conser- 
vative force. As we shall see below, the conservative force must 
be included in the calculation when the evolution of the orbital 
phase is required to stay accurate to order e'' during a radiation- 
reaction time. In a context where the radiation-reaction force 
would contain a second-order term, the same accuracy would be 
achieved only after including this term as well in the calculation. 
We thank Tanja Hinderer and Eanna Flanagan for making this 
point clear to us. 
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C. Multi-scale approximation: time 

We now wish to construct a multi-scale approxima- 
tion for the time function t(0). To begin we recall that 
Eq. p.l2p must be integrated numerically even in the 
unperturbed situation, when p, e, and uj are all constant; 
the approximation, therefore, will also involve a numer- 
ical integration. A lazy option presents itself: The time 
function could be obtained simply by inserting our multi- 
scale approximations for p (</>), e(0), and i-u{4>) within 
Eq. ()3.12|) and performing the integration numerically. 
In an effort to obtain maximum analytical insight, how- 
ever, we choose to proceed differently. 

We substitute Eqs. (|3.13p - (|3.15p and the explicit ex- 
pressions of Eqs. (jX^ - f?:^ into Eq. (jXT^ . and we 
expand in powers of e. Through first order we obtain 

3/2 

./ _ Po 



(1 -I- eo cosu)^ 

1 1/2 4 -t- Co cos u 



2 ''''''^ (1 -I- ep cosw)^ 

1 sinu -|- eo sin2u 

2 (1 -I- eocosu)'' 



(3.32) 



where v = 4> — ujq, and po, eo, and loq are the functions of 
z := Cr r^ di splayed in Eqs. (j3?29l) - ((3?3T|) . The first hnc 
in Eq. p.32p is obtained by inserting p = po, e = eo, and 
UJ = ujQ within Eq. (j3.f 2^ : the second and third lines are 
contributed by the oscillatory terms pi, ei, and wi. 

We wish to find an approximate solution to Eq. (|3.32p . 
and once more we wish to distinguish between secular and 
nonsecular terms. Two sources of complications present 
themselves. First, while it was easy in Eq. (|3.16[) - (|3.18[) 
to separate the secular terms from the oscillations, the 
factors of 1 -f- eo cos v in the denominators of Eq. (|3.32p 
make this separation more difficult. Second, while p, e, 
and UJ, are simply constant at the unperturbed level, the 
Keplerian version of t{(l)) is already a complicated func- 
tion of (j) that contains secular and oscillating terms. The 
situation here is therefore more complicated, but we will, 
nevertheless, be able to express the solution to Eq. (|3.32p 
in the form 



t(0) =to(0) +eti(0) 



(3.33) 



with to((/>) incorporating the Keplerian behavior (includ- 
ing secular terms and oscillations) as well as the secular 
changes produced by the electromagnetic self-force, and 
with eti((/)) being purely oscillatory. 

To isolate the oscillatory terms in Eq. (|3.32p we calcu- 
late the averages 



2,71 



J{4>')d4>' 



0— 7r 



of the various functions of (j) that appear on its right-hand 
side; these averages are calculated while keeping po, eo, 
and ujQ constant over the integration domain. Defining 



/i = (1 + eo cosu) 2, /2 = (4 + eocos'i;)(l-|- eocosi;) ^, 
and /s = (sinu + eosin2w)(f + eocosw)"^, we find that 
= (l~e2)-3/2, = i(8 + e2)(l-e2)-5/2, and 

(/a) — 0. The function multiplying in Eq. (|3.32p is 
therefore purely oscillatory, but the function multiply- 
ing Ec contains a secular component. To isolate this we 
rewrite Eq. (I3.32p into the equivalent form 



t' - 



3/2 

Po 



(I -f eo cosw)^ 



1 Po 

4 ' 
1 



4po(l-e2) 
y^eo 9eo -I- 3(2 -f e^) cos v 



sm V 



r^rreo- 



(I + eo cosi;)3 
eo sin 2v 



I -I- eo cosvY 



(3.34) 



in which a term jtcPxi^i^ + e§)(I — eo)^"^(l -f eo cosu)~^ 
was removed from the second line in Eq. (|3.32p and in- 
serted within the first line. In Eq. (|3.34p . the functions 
that appear in the second and third lines are purely os- 
cillatory. 

It is important to notice that the term proportional to 
Ec in the first line of Eq. (|3.34p is a secular correction to 
t' that originates with the oscillatory terms pi, ei, and 
uji in the orbital elements. These oscillations combine 
in a nonlinear fashion, and they contribute an additional 
secular term beyond the one that comes from po, eo, and 
ujQ. It would be a significant mistake to discard the os- 
cillations in the orbital elements when constructing the 
time function. 

The solution to Eq. ([glM)) is 



3/2 

Po 



(1 + eoCOSw')^ 



'4po(l-e^) 



1 pJ/^eo9eo + 3(2-)-e^)cosv' 

4'^ 



1 



-0 



(1 -I- eo cosv'Y 



1 



eo- 



sinu' -t- eo sin 2v' 
(1 + eo cosw')'^ 



(3.35) 



where v' :— 



ujo{(p'). We must leave the first integral 
alone, but we shall manage to evaluate the second and 
third integrals. Because the changes in po, eo, and ujq 
are of order e, because the second and third integrals 
already come with a factor of e in front, because the 
integrands are purely oscillatory functions, and because 
the calculation of t{(j)) is carried out consistently to first 
order in e, we are permitted to treat po, Cq, and ujq as 
constants when evaluating the integrals. We thus obtain 



3/2 
Po 



(1 + eo cos v')^ 



'4po(l-e2) 



1 /'^ 

1 pp'-eo 3 sinu + |eo sin2w 

2*^''l-e2 Cl -U^„r.nc,A2 



1 



(1 + eo cosi;j 
4eo cos V 3 



■4en 



(l + eocosu)2 {1 + eoy 



(3.36) 
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We have not yet achieved the form of Eq. (|3.33p . The 
reason is that while the function in the second line of 
Eq. (|3.36p is purely oscillatory (it has a zero average), this 
is not true of the function in the third line. Defining — 
(3 + 4eocosw)(l + ecosw)"^, we find that (/4)0 = (3 — 
4eo)(l — Cq)^"^/^ and this combines with the second term 
on the third line to contribute secular terms. Removing 
these from the third line and inserting them within the 
first line, we finally arrive at the desired expression for 
the time function t (</)). 

Our final result is that t((/)) can be expressed as in 
Eq. (0331), with 



to 



3/2 

Po 



Q 1^1 + Co cos u')^ 
3 + 4eo 



eg 



■4po(l 



1 



4e 



2 1 



(l + eo)2 (l-e2)3/2 



(3.37) 



and 



eti 



1 /2 S 

1 pg' eo 3sinu + ^eo sin2z; 
'2'^'=l-e2 (1 + eocosi;')2 
1 r 3 + 4eo cos V 3 — ie^ 



(l + eocosi;)2 (1 - 6^)3/2 



(3.38) 



where v = (f> ~ luq, and po , eo , and luq are the functions 
of z := e„(j) displayed in Eqs. (|3.29p - (|3.3ip . By design, 
to(0) incorporates the Keplerian behavior (including all 
Keplerian oscillations) in addition to the secular changes 
produced by the perturbing force; the function eti(0) is 
purely oscillatory, in the sense that its 0-average is zero. 
A comparison between the exact time function t{(j)) and 
the multi-scale approximation is presented in Fig. 5. 



IV. DISCUSSION 

With the technical details out of the way, we may now 
return to the themes introduced in Sec. I. To launch the 
discussion we summarize the main results obtained in 
Sec. III. 



A. Summary of our results 

In the method of osculating orbital elements, the mo- 
tion of a charged particle subjected to the electromag- 
netic self-force of Eq. (|1.2p is at all times described by 



K0) = 



p 



1 



(4.1) 



the Keplerian relation of Eq. (|3.ip . The orbital elements 
p, e, to, however, acquire a (?!)-dependence that accounts 
for the perturbation created by the self-force. With the 
definitions of Eq. (|3.8p , these quantities evolve according 
to Eqs. (|5?g)) - ([5TT|l . and integrating Eq. ((XT^ produces 
t{4>). The motion is then fully determined. 




20 30 
(p / (2 It) 



FIG. 5: Multi-scale approximation for t{<j>/2TT) compared with 
exact numerical results. The numerical conditions are the 
same as in Fig. 1. The solid curve in red shows the exact 
evolution as computed numerically. The dotted curve in blue 
shows the evolution as predicted by the multi-scale approxi- 
mation to + eti, where to incorporates the Keplerian behavior 
in addition to the secular changes produced by the perturb- 
ing force, and where et-i{(j>) is purely oscillatory. The dashed 
curve in green is a plot of to only. The large panel shows the 
entire evolution from (j> = to <j> — lOOvr. The first inset (top 
left) shows the evolution in the small interval 6 < 0/(2^) < 9. 
The second inset (bottom right) shows the evolution in the 
interval 41 < (I>/{2tt) < 44. 



The evolution equations can be integrated numerically, 
or they can be integrated analytically via a multi-scale 
analysis that produces a faithful approximation over the 
long interval < < e~^. Moreover, the multi-scale 
analysis produces a clean separation of the solutions into 
secular and oscillatory pieces. The secular changes in the 
orbital elements are described by the functions po, cq, and 
ujQ displayed in Eqs. p.29p - (|3.3ip . We copy them here 
for convenience: 



Po 
eo 



P*(l-3e„ 

4 + e 



2/3 



l + Cc 



4e* 

8 + 5e* + 4e*2 



1 - 3e,.i.(; 



.1/2 



2e„ 



4e* 
1 - (1 - 3er 



,1/3 



where 



err 




(4.2) 
(4.3) 

(4.4) 

(4.5) 
(4.6) 



and where p* := p{(j} = 0) and e* :— e{(j3 = 0); we have 
set cl'((/) = 0) = 0. The oscillatory changes in the or- 
bital elements are given by epi, eei, and ecui displayed in 
Eqs. (IS2ni)-(E2Sl)- 
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The piece of the time function that incorporates Ke- 
plerian behavior and secular changes produced by the 
electromagnetic self- force is and this is displayed in 
Eq. (|3.37p . We copy it here for convenience: 



^0 = \^ 




Co cosw')^ 



1 



(8 + eg)p* 
^4po(l-e2) 



3 + 4eo 



3-4eg 



(l + eo)2 (l-e2)3/2 



(4.7) 



where 



ujo{(f>')', this, like the Keplerian time 



function, is expressed in terms of an integral that must 
be evaluated numerically. The oscillatory piece of the 
time function is eti, and this is given by Eq. p.38p . We 
recall that the term proportional to Cc inside the integral 
is produced by the oscillatory pieces pi, ei, and uji of the 
orbital elements; the oscillations combine to produce a 
secular correction to the time function. 



B. Conservative and dissipative terms in the 
self-force 



The effect of each term in Eq. (|1.2p can easily be iden- 
tified if we focus our attention on the secular changes 
in the orbital elements and time function described by 
Eqs. (|4.2p - (|4.7p . Because these accumulate in the long 
run, while the oscillations that are not contained in those 
equations average to zero, it is clear that it is the secular 
pieces of e(0), w(0), and t{(p) that arc the most 

important to capture. 

The effects of the conservative piece of the self-force 
are identified by selecting the terms in e^, the effects of 
the radiation-reaction piece are identified by An ex- 
amination of Eqs. (|4.2p - (|4.7p reveals that the radiation- 
reaction force drives secular changes in the principal or- 
bital elements p and e, but that it affects w only indi- 
rectly, and only if there is a conservative force. In ad- 
dition, the radiation-reaction force affects the time func- 
tion indirectly through the changes in the principal ele- 
ments, and also directly as can be seen in the second line 
of Eq. (|4.7p . On the other hand, the conservative force 
drives secular changes in the positional element lu, but 
it affects e only through the factor that comes in front 
of (1 — 3err</')^^^- In addition, the conservative force af- 
fects the time function directly, as can be seen from the 
correction term inside the integral. 

The combined effects of the conservative and radiation- 
reaction forces on the time function can be neatly sum- 
marized by computing P := J^^'^'^ t{4>') d(j)' , the period 
of an orbital cycle. Ignoring the changes in the orbital 
elements while performing the integration, we obtain 



P = 2ti^ 



M{1 



1 



1 p^S_ 

4'%T 



+ e 



(4.8) 



The factor in front of the large brackets is the Keple- 
rian period expressed in terms of the changing orbital 



elements; these changes, we recall, are driven by the 
radiation-reaction force. The second term gives the cor- 
rection contributed by the conservative force. It may be 
recalled that this correction originates from oscillations 
in p, e, and and it may be noted that it becomes large 
when eo — >■ 1. 



C. Limitations of the radiative approximation 

As we have defined it in Sec. I, the radiative approxi- 
mation is obtained by setting ec = in our results. This 
preserves the secular changes in p and e, but it com- 
pletely turns off the secular evolution of w. In addition, 
the radiative approximation discards an important cor- 
rection term in the time function, the one proportional 
to ec in Eq. ()4.7p . This term, in fact, dominates over 
the radiation-reaction corrections, because Cc is numeri- 
cally larger than Crr by a factor of order ^p* jM ^ 1, as 
can be seen from Eqs. (|4.5p and (|4.6p . In addition, we 
have noted that the correction becomes increasingly large 
as eo increases toward unity. As a result, the radiative 
approximation provides a poor estimation of the orbital 
period, as can be seen from Eq. ()4.8p . 

The combined effect of omitting the secular changes in 
Lo and missing an important correction in the time func- 
tion can be seen most clearly by examining the radial 
phase variable $ :— </> — w expressed as a function of time; 
this is the function that appears in r(t) = p/(l + e cos $), 
as can be seen from Eq. (|4.ip . In Fig. 6 we compare the 
results of two numerical computations, one carried out 
with the full electromagnetic self-force (including con- 
servative and radiation-reaction pieces), the other car- 
ried out in the radiative approximation, with only the 
radiation-reaction piece of the self-force. The figure, and 
the quantitative analysis presented in the caption, re- 
veal very clearly that the radiative approximation gives 
a rather poor representation of the phase function; in our 
simulation, the mismatch after 50 orbits is nearly three 
full radial cycles. The origin of the discrepancy is also 
clearly identified in the caption: It is the missing conser- 
vative correction to the time function ^(0) that is mostly 
responsible for the phase mismatch. 



D. Secular approximation: ^-average 

The multi-scale approximation method of Sec. Ill was 
adopted precisely because it produces a clean separation 
between secular and nonsecular terms in the expressions 
for the orbital elements and the time function. The 
zeroth-order quantities displayed earlier in this section 
were constructed so as to represent the secular changes, 
and the oscillatory corrections were carefully designed to 
average to zero. It is clear, therefore, that Eqs. (|4.2p - 
(|4.7p achieve the goals of a secular approximation, and 
we would be justified to write 



Pscc(0) =P0 (</>), 



:(0) = eo((?!)). 
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FIG. 6: Effect of the radiative approximation on the radial 
phase function where $ = cj> ~ uo{(j)). Plotted are <l?/(27r) 
versus t = ^ M /p*^t, with the same numerical conditions as 
in Fig. 1. The lower curve in red is the exact evolution as com- 
puted numerically with the help of the full self-force, which in- 
cludes conservative and radiation-reaction pieces. The upper 
curve in blue is the evolution obtained in the radiative approx- 
imation, in which the conservative force is switched off. At the 
end of the integration, for t = 875, we have <l>/(27r) = 50.40 
under the action of the full self-force, with a contribution 
uj/{2-k) — —0.073 coming from the periapsis advance. The 
radiative approximation gives instead $/(27r) — 53.06, with 
a contribution uj/{2'k) = 0.00016 coming from the shift in pe- 
riapsis. The total phase mismatch is A$ = 2.66(27r), nearly 
three radial cycles out of fifty orbits. Because the contribu- 
tion from uj is small in both cases, we conclude that the main 
source of error is contained in the missing conservative terms 
in *(<?!>). 



(4.9) 



This, in the language introduced in Sec. II, consists of 
defining the secular approximation by averaging the ex- 
act solution (as represented by the multi-scale approxi- 
mation, which has been demonstrated to be faithful to 
the exact numerical results) over the orbital parameter 

The question we wish to explore here is whether the 
secular approximation of Eq. (|4.9p could be formulated 
directly, without the help of the exact solution. The dy- 
namical equations that govern the secular approximation 
can be obtained by differentiating Eqs. (|4.2p - (|4.7p with 
respect to (j). We obtain 



Pscc 

e' 

t' 



*3/2n 



(errP"'~jesocPsoc 
(ecP*)ftcc: 



3/2 



VpIoc/m 



(1 



; COS vy 



1 



(eeP*)(8 + eL) 



4(1 



(4.10) 
(4.11) 

(4.12) 
(4.13) 



in which v = 4> — a;soc(0) and where we have ignored (as 
we should) terms of order . These equations must come 



Pscc \ 
^scc ( 

(^xnr ( 



= 0) = p*, 
^0) ^ 



0) 
0) 



1 + e, 



4 -He* 
' 4e* 
5e* + 4e*2 

4e* ' 



3 + 4e* 
(l + e*)2 

3 - 4e*2 
" (1 - e*2)3/2 



(4.14) 
(4.15) 

(4.16) 



(4.17) 



in order to reproduce precisely the secular evolution pre- 
dicted by the multi-scale approximation. 

The differential equations for Pscc, escc, and Wsec are 
easy to motivate: Eqs. (|4.10[) - (|4.12[) are the same as 
Eq. ()3.19|) . and they can be obtained directly by sub- 
mitting the exact equations (I3.9p - ()3.1ip to an averaging 
procedure. The differential equation for igcc, however, is 
not so easy to justify. It is not reproduced by averaging 
Eq. (|3.12p over </>, which would fail to account for the 
important conservative correction proportional to e^.; the 
averaging would also remove the Keplerian oscillations of 
the time function. And the initial values of Eqs. (|4.14p - 
(|4.17p cannot be justified at all without knowledge of the 
oscillatory terms in the multi-scale approximation. To 
integrate the differential equations with the approximate 
initial conditions Pscc(O) = p* , escc(O) = e*, Wscc(O) = 0, 
and tsec{0) = would produce solutions that are offset 
from the exact solutions by quantities of order e. It is 
noteworthy that while the corrections to Pscc(O), Wscc(O), 
and iscc(O) come as additive terms that become increas- 
ingly irrelevant as increases, the correction to esoc(O) 
comes as a multiplicative factor, this correction never be- 
comes irrelevant. 

Our main conclusion is this: The secular approxima- 
tion defined by the system of equations (|4.10p - (|4.17p 
would be very difficult to formulate without prior knowl- 
edge of the exact solution, as represented by the faithful 
multi-scale approximation. This conclusion reflects the 
lesson learned from the illustrative example described in 
Sec. II. 



E. Secular approximation: t-average 

The secular approximation considered in the preceding 
subsection is obtained by removing the oscillations in 
from the exact expressions for the orbital elements. Be- 
cause we are ultimately interested in the time behavior 
of the elements, it is perhaps more meaningful to define 
the secular approximation by averaging with respect to 
t instead of 0. In this alternative secular approximation, 
we write 



Pscc 



:c(0) = (a;)t(0). 



icc(0) 



(e>t(0), 

- {t)M) 



(4.18) 
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in place of Eq. (|4.9p . where the time average is defined 
as in Eq. 



j';^:qmdt/d^')d^' 



st^:{dt/dcp')d^' 



(4.19) 



The oscillations contained in dt/d(t), as revealed in 
Eq. p.32p . ensure that this version of the secular ap- 
proximation is quite distinct from the version examined 
in the preceding subsection. 

Performing the calculations produces 



Pscc 
^scc 

^sec 



4 eopo 



LOq ■ 



hel ~ 2(1 



elf' 



2 3/2 



(4.20) 
(4.21) 

(4.22) 



It is clear that these expressions do not agree with those 
of the preceding subsection, in which we made the as- 
signments Psec = PO, Gsoc = eo, and Wscc = Wq. 

In view of the complexity involved, we shall not at- 
tempt to find an explicit expression for tsec ■ Nor shall we 
prolong the discussion by writing down dynamical equa- 
tions and initial conditions for Pscc, escc, ^sca and tgoc- 
We can simply jump to the main conclusion, which is 
the same as in the preceding subsection: The dynami- 
cal equations and initial conditions associated with this 
version of the secular approximation would be very dif- 
ficult to formulate without prior knowledge of the exact 
solution. Once more this conclusion reflects the lesson 
learned from the illustrative example described in Sec. II. 



V. CONCLUSIONS 

We examined the motion of a charged particle in a 
weak gravitational field. In addition to the Newtonian 
gravity g exerted by a large body of mass A/, the particle 
is subjected to the electromagnetic self-force described by 
Eq. p.2|) . As we have argued in Sec. I, this toy problem 
shares many of the features of the gravitational self-force 
problem, and yet it is sufficiently simple that it can be 
solved completely with simple numerical methods, and 
virtually completely with simple analytical methods. 

After subjecting the equations of motion to a multi- 
scale analysis in Sec. Ill, we summarized our main results 
in Sec. IV and investigated the main themes introduced 
in Sec. I. We first examined the roles of the conserva- 
tive and radiation-reaction pieces of the self-force. We 
showed that the radiation-reaction force drives secular 
changes in the principal orbital elements p and e, while 
the conservative force drives secular changes in the po- 
sitional element lj as well as in the time function t(0). 
This led us to our first conclusion: 

The radiative approximation to the true self- 
force does not account for the secular 



changes in all the orbital elements; this 
gives rise to an important phase mismatch 
between an orbital evolution driven by the 
radiation-reaction force, and one driven 
by the true self-force. The radiative ap- 
proximation does not achieve the goals of 
a secular approximation. 

This was also the conclusion of our previous work (pa- 
per I: Ref. [l3|), but we believe that we have established 
these statements more firmly in this work. In addition, 
the source of the phase mismatch was correctly identified 
here, while it was attributed incorrectly in paper I: it is 
the conservative correction in the time function that is 
mostly responsible for the dephasing, and not the secular 
change in uj. 

We next considered the issue of formulating secular 
approximations to the dynamical equations that govern 
the evolution of the orbital elements. Having access to 
a faithful, analytical representation of this evolution, as 
provided by the multi-scale approximation, it was an easy 
task to perform averages and to obtain expressions that 
capture the secular changes in the orbital elements (and 
the time function). And having access to those expres- 
sions, it was again an easy task to identify the differ- 
ential equations that govern their behavior, as well as 
the appropriate initial conditions. The issue, of course, 
is whether the simplified dynamical equations, those that 
would govern the purely secular changes in the orbital el- 
ements, could be obtained directly in a context in which 
the exact solutions are not known. Our answer is in the 
negative, and this led us to our second conclusion: 

A secular approximation to the exact differ- 
ential equations and initial conditions, de- 
signed to capture the secular changes in 
the orbital elements and to discard the 
oscillations, would be very difficult to for- 
mulate without prior knowledge of the ex- 
act solution. While some of the approx- 
imate differential equations can be ob- 
tained by submitting the exact equations 
to an averaging procedure, other equa- 
tions cannot be obtained so simply. And 
even if the correct differential equations 
can be identified, their integration must 
proceed from initial conditions that differ 
from the exact initial conditions; the dif- 
ference is determined by the oscillations, 
and those must be known before the ap- 
proximate initial conditions can be pre- 
scribed. 

In addition to these issues, the formulation of a secu- 
lar approximation must resolve a fundamental ambiguity: 
Which oscillations are to be removed? In our analysis we 
had to distinguish carefully between taking a 0-average to 
remove oscillations in and taking a ^-average to remove 
oscillations in t. Different choices lead to different secu- 
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lar approximations, different dynamical equations, and a 
different prescription for initial conditions. 

The work presented here leaves a number of questions 
to be examined. The most important one is this: Do 
the conclusions of this paper have any relevance to the 
gravitational self-force? While our analysis of the elec- 
tromagnetic self-force leaves no room for controversy, the 
question of how our results will transfer to the more inter- 
esting case of the gravitational self-force might be cause 
for debate. We believe that the analogy between the 
electromagnetic and gravitational self-forces is close, we 
believe that our general conclusions do carry over to this 
case, and we believe that our work serves as a useful 
cautionary tale for the gravitational self-force. But we 
admit that the analogy relies on the usual formulation of 
the gravitational self-force in the Lorcnz gauge, and that 
the analogy may be lost in alternative formulations — 
the gravitational self- force is not gauge invariant, and its 
effect on the description of orbital evolutions will dep end 
on the choice of gauge. For example, Mino 0, l22l | has 
proposed a formulation of the gravitational self-force in 
a "radiation-reaction gauge" in which the full self-force is 
equal (for a radiation-reaction time) to the radiative self- 
force. In Mino's proposed formulation, the radiative ap- 
proximation is exact over a radiation-reaction time, and 
the issues raised here may not at all be relevant. How 
our conclusions might apply to the gravitational case is 
indeed a controversial topic, but we consider its discus- 
sion to be beyond the scope of this work. Indeed, this 
paper is concerned with the electromagnetic self-force, 
and the case of the gravitational self-force is considered 
separately in a companion paper (23j . In our companion 
work we argue that the Lorenz-gauge formulation of the 
gravitational self-force is physically meaningful, that the 
Lorenz gauge is most likely to keep quantities other than 
the self-force (such as the gravitational potentials) under 
control, and that the conclusions of this paper do carry 
over to the gravitational case. 
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APPENDIX A: KEPLERIAN MOTION 

In this Appendix we provide a complete description of 
Kepler's problem. This material is well known, and can 
be found in any textbook on celestial mechanics (see, for 
example, Ref. 24]), but we include it here for complete- 
ness and as a way of defining our notation. 

Two bodies of masses mi and TO2 move under their 
mutual gravitational attraction. The equation of motion 
for the relative position r := ri — r2 is 



where a := cPr/dt^ is the relative acceleration vector, 
and g — —Mr/r^ is the gravitational field. Here M = 
nil + TO2 is the total mass, and r — \r\ is the distance 
between the two bodies. We set G — 1. 

Conservation of angular momentum implies that the 
motion takes place within a fixed plane. We use polar 
coordinates (r, (/)) in this plane, and we resolve all vec- 
tors in the associated basis (r,<^). The relation with 
the Cartesian description is a; = rcoscf), y = rs'mcj), 
f = cos ipse + s'm(py, and cj) — —simpx + cos (f>y. 
The position vector is r — rf, the velocity vector is 
V = f r + r(j) 4>, and the acceleration vector is 



a = (f - r0^)r -I- -^{r^(l))4>. 

r at 



1 d 



(A2) 



An overdot indicates differentiation with respect to t. 
Equations (|XT|) and imply 



r'^(p — constant \/Mp, (A3) 
which defines the semilatus rectum p. We also have 
.. M Mp 



= 0, 



which integrates to 



1.2 

r 



M 
r 



Mp 
2^ 



constant 



(A4) 



(A5) 



The constant is the system's conserved energy per unit 
reduced-mass, and the last equation defines the eccen- 
tricity e. 

Eliminating time from Eqs. (|A3[) and (jA4[) produces a 
differential equation for r{(j)) which integrates to 



P 



I + e cos(0 — Lu) 



(A6) 



where lo is an additional constant of the motion. This 
equation describes an off-centered ellipse of semi-major 
axis 



P 



1 



(A7) 



and eccentricity e. The constant uj, known as longitude of 
periapsis, determines the orientation of the ellipse in the 
plane. The orbit is at periapsis r = p/{l + e) whenever 
cos((/) — oj) ~ 1 , and is at apoapsis r = p/(l — e) whenever 
cos((/) — u) ~ —1. 

Equations (jA3|) and (jA6| imply 



(A8) 




and 



9, 



(Al) 



■ [l + ecos(0 — uj)]' 



(A9) 
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This last equation integrates to 




d(t)' 



[l + ecos(0' — Lu)]' 



(AlO) 



and determines the time. The fourth (and final) constant 
of integration tpon is time at periapsis, and is such that 

t{(f) = uj) = t 

period is 



pen- According to Eq. (jAlOp the orbital 

(All) 



n 



where n is known as the mean motion. 



APPENDIX B: OSCULATING ORBITAL 
ELEMENTS 

In this Appendix we develop a method of osculating or- 
bital elements for the integration of the equations of mo- 
tion associated with a perturbed Keplerian orbit. The 
general idea is very old, and many variations of this 
method can be found in the literature (see, for exam- 
ple, Ref. [13|)- But we find that the version presented 
here is perhaps a little unusual, while being especially 
convenient and well suited to our purposes. For these 
reasons we judge it worthwhile to develop it in full here. 

We consider the equations of motion 



o. = g + f, 



(Bl) 



in which / is a perturbing force (divided by the system's 
reduced mass) that depends on the relative position vec- 
tor r and (possibly) the relative velocity vector v. (The 
notation is introduced in Appendix A.) We seek to in- 
tegrate Eq. (|B1[) for r{t) using a method of osculating 
orbital elements. We assume, for simplicity, that the per- 
turbing force can be decomposed as 



f = Bf + Scl), 



(B2) 



so that it lies within the orbital plane. The perturbed 
orbit, therefore, will stay within the same plane. 



Let 



First formulation 



:= {p, e,cj,tpori} 



(B3) 



collectively stand for the Keplerian orbital elements in- 
troduced in Appendix A, let 

rK{I^,t) (B4) 

stand for the position vector of a Keplerian orbit, and let 

VKiI^,t) (B5) 



be the Keplerian velocity vector. The method of oscu- 
lating elements states that the perturbed motion is de- 
scribed at all times by Eqs. (|B4[) and (|B5[) . but that the 
orbital elements acquire a time dependence. In mathe- 
matical terms, the position vector of the perturbed orbit 



r^rK{l^{t),t) (B6) 
and its velocity vector is 

v = VK{l'^{t),t). (B7) 
Differentiating Eq. (|B6p with respect to time yields 



dvK dl^ drK 
" dP^~dr ^ ~W 

The second term, in which tk is differentiated while keep- 
ing constant, is recognized as VK^ the Keplerian ve- 
locity vector. Comparing with Eq. jBTj gives 



^1^ = 



(B8) 

Differentiating Eq. (jB7p with respect to time yields 

i9i)K dl^ dvK 
dl^ dt '^~dr' 

The second term gives g, and comparing with Eq. (jBip 
gives 



dvK 



(B9) 



Equations (jBSp and (jB9|l can be solved for in terms 
of the perturbing force. The equations of motion have 
become a system of first-order differential equations for 
the orbital elements. The method of osculating orbital 
elements therefore transforms the original phase space 
spanned by (r, v) into a new phase space spanned by the 
coordinates I^. In the planar context considered here, 
the original phase space is spanned by (r, (jj, f, 0) while 
the new phase space is spanned by (p, e, lu, ipcri)- 
Concretely the equations of motion are 



M 



dt 



rS. 



(BIO) 



By virtue of Eq. (jA3p and the osculating conditions of 
Eqs. jBGl and jBTl, r^(j) = ^/Mp and the second of 
Eqs. (|BT0)) implies rS ^ \^fMjpp. Inserting Eq. 
yields 



1 



M 1 + e cos(0 — lu) 
the new equation of motion for p{t) . 



(Bll) 
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To work out the remaining equations we substitute 
Eq. jM]) into the first of Eqs. (|BT0)) . This gives 



where 



e M , 





(B12) 



after cancehng out all Keplerian terms. An additional 
equation is obtained by differentiating Eq. (|A6|) with re- 
spect to time and demanding that the result be compat- 
ible with Eq. (jASp . After some algebra we obtain 



pcos( 



- w 



ep sin( 



= p 

1 + e cos((/) — cj) 1 + e cos((/) — w) 
Equations (|BTT|) - (|BI3)) imply 



uj. (B13) 



sm 



V_ 

M 

e + 2 cos( 



uj)R 

' — uj) + e cos^(0 — uj) 



1 + ecos((/) — Lu) 



(B14) 



and 



cos 



sin((/) — cj) [2 -I- e cos((/) — ijj)] 
1 + e cos{(j) — uj) 



S 



(B15) 



In these equations, </) is a function of time that must be 
obtained by integrating Eq. (jA9p . 




ecos(0 — u)]' 



in which p, e, and w are now time-varying orbital ele- 
ments. 

Our system of equations currently leaves out tperi, the 
fourth orbital element. An equation for iperi, however, 
will not be required. 



2. Second formulation 

The preceding system of equations achieves a cleaner 
structure if we change the independent variable from t to 
(j) via Eq. (|A9|) . Writing, for example, p' := dp/d(p — p/<j), 
we obtain 



P 



2p^ 



1 



M (1 + ec)3 

p'' 

M 
M 



(1 + ec)2 
c 



(1 + ec)2 
1 

M (1 + ec)2 ' 



5", 

e + 2c + ec2 ' 
^+ (l + ec)3 



(l + ec)3 



(B16) 
(B17) 
(B18) 
(B19) 



c := cos{(f> — uj), s :— sin((/) — u). 



(B20) 



The first three equations for e((/)), and consti- 
tute a closed system that can be solved independently of 
the fourth equation, which determines These equa- 
tions are exact, they are convenient to deal with, and they 
can easily be implemented numerically. (The equations 
are ill-behaved when e — > 0; a transformation to new 
variables a = ecosw, f3 = esinw eliminates this pathol- 
ogy.) It is understood that the system of equations is 
accompanied by the Keplerian representation of the mo- 
tion, that is, equations such as r(0) = p/[l-t-ecos(0 — w)] 
and r' = epsin(0 — uj)/[1 + ecos(0 — uj)]'^. 

The second formulation of the method can be under- 
stood as follows. Let 



■.^{p,e,uj} 



(B21) 



collectively stand for the relevant orbital elements, let 

rK(/-^,0), iK(/-^,0) (B22) 

stand for the position vector of a Keplerian orbit, param- 
eterized by longitude 0, and let 



del) ' 



(B23) 



The Keplerian velocity vector can then be expressed as 

The method of osculating elements states that the per- 
turbed motion continues to be described by Eqs. (jB22[) 
and (|B23p . but that the orbital elements acquire a 0- 
dependence. In mathematical terms, the position vector 
of the perturbed orbit is 



and we impose also 



The first two equations are equivalent to Eq. (jB6p . and 
the last two equations are equivalent to Eq. (|B7p . The 
second of Eqs. (|B25p is the same as Eq. (IB19|) . 
Differentiating Eq. (|B24p with respect to (f> yields 



t = t^{l^{4>),<f) (B24) 



t' = t'^{l^{^),<t>). (B25) 



drK dl^ 



r = 



dl^ d(j) ^ d(f)' 
Comparing with Eq. (IB25|) gives 



Differentiating Eq. (jB7p with respect to 
by t' from Eq. (|B%1) yields 



(B26) 
and dividing 



V 



dvK dl^ dvu 



dl^ d4> d4) 
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The second term gives g, and comparing with Eq. (|B1 
gives 



1 dv 



T^A - J- 



(B27) 



Equations (IB26|) and (jB27|) can be solved for in terms 
of the perturbing force, and the end result is the system 
of Eqs. (jB16jl - (|B18|l . In this formulation the method of 
osculating orbital elements transforms the original phase 
space spanned by r{t), (l>{t), f{t), and <j){t) into a new 



phase space spanned by p{4>), e(0), w ((/>), and t{(p). 



The second formulation of the method of osculating 
elements is distinguished by the facts that it involves (p as 
a running orbital parameter, and it removes iperi from the 
list of phase-space variables. This formulation leads to 
the important advantages that Eqs. (|B16|) - (|B18|1 form a 
closed set of equations; these equations can be integrated 
first, and t{(p) can be recovered at a later stage by solving 
Eq. (IBT91) . 
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